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Using renormalization group techniques we have derived analytic formulae for the next-to-leading order heavy- 
quark coefficient functions in deep inelastic lepton hadron scattering. These formulae are only valid in the 
kinematic regime ^ m^, where and stand for the masses squared of the virtual photon and heavy 
quark respectively. Some of the applications of these asymptotic formulae will be discussed. 



1. INTRODUCTION 

Deep inelastic electroproduction of heavy 
flavours is given by 



-(£i)+P(p) 



"(^2)+Q(pi)(Q(p2))+'X'.(l) 



When the virtuality -q^ ^ Q'^ > { q = ii - £2 ) 
of the exchanged vector bosons is not too large 
( ^ Mf ) the above reaction only gets a 
contribution from the virtual photon and we can 
neglect any weak effects caused by the exchange 
of the Z-boson. If the process is inclusive with 
respect to the hadronic state 'X' as well as the 
heavy flavours Q{Q), the unpolarized cross sec- 
tion is given by 



dxdy 
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-y'^^L(x,Q^m2)l , (2) 
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where S denotes the square of the cm. energy 
of the electron proton system and the variables x 
and y are defined as 

x = ^{^<x<\),y^^- (0<y<l),(3) 



with 

- = = xyS 



(4) 



In the QCD improved parton model the heavy 
flavour contribution to the hadronic structure 
functions, denoted by Fi{x,Q'^ ,171^) {i = 2,L), 
where m stands for the heavy quark mass, can be 
expressed as integrals over the partonic scaling 
variable. This yields the following results 



Fi{x,Q'^,m'^) =x 
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where z — Q"^ /{s + Q^) and s is the square of the 
photon-parton centre-of-mass energy (s > Am?). 
Here the upper boundary of the integration is 
given by z^ax = Q'^/{Am'^ + Q^). 

The function G{z, jj?') stands for the gluon den- 
sity whereas the flavour singlet and flavour non- 
singlet combinations of the quark densities are 
given by T,{z,ix^) and A{z,fi'^) respectively. In 
the above expressions the charges of the light 
quark and the heavy quark are denoted by and 
en respectively. Furthermore, n/ stands for the 
number of light quarks and n denotes the mass 
factorization scale, which we choose to be equal 
to the renormalization scale. The latter shows 
up in the running coupling constant defined by 
Usi^J?)- The heavy quark coefficient functions are 
denoted by and Hi j (i = 2, L; j = g). The 
distinction between them can be traced back to 
the different photon-parton production processes 
from which they originate. The functions Lij 
are attributed to the reactions where the virtual 
photon couples to the light quark, whereas the 
Hi j originate from the reactions where the vir- 
tual photon couples to the heavy quark. Hence 
Lij and Hij in (5) are multiplied by ef and eg 
respectively. The superscripts NS and S on the 
heavy quark coefficient functions refer to flavour 
non-singlet and flavour singlet respectively. Fur- 
thermore the singlet quark coefficient functions 
if ,j and Hf^^ can be split into non-singlet and 
purely singlet (PS) parts, i.e.. 



-NS 



-PS 



H. 



NS 



PS 



(6) 



(7) 



with -fffg = in all orders of perturbation theory. 

In ||l[ the heavy quark coefficient functions Lij 
and Hi j have been exactly calculated up to a^. 
Expanding them in a power series in {as/ATT)'^ 



they receive contributions from the following par- 
ton subprocesses 

l*iq)+9iki)~^Q{pi) + QiP2), (8) 
7*('Z) + gih) ^ 9{k2) + Q{pi) + Q{P2) , (9) 
and 

l*{q) + qmki) ^ qmk2) + Q{pi) + Q{P2) .(10) 

For reaction (9) one has to include the virtual 
gluon corrections to the Born process (8). The 
contributions from (8) and (9) to the heavy quark 
coefficient functions are denoted by h[]^ and 
respectively. The parton subprocess (10) has two 
different production mechanisms. The first one is 
given by the Bethe-Heitler process (see figs. 5a, 5b 
in [Ql) leading to Hf^'^^^ and the second one can 
be attributed to the Compton reaction (see figs. 
5c, 5d in ||l|). Notice that Lf^^ and Lf^ are zero 
through order a^. Then, from (6), one infers that 
^NS,(2) _ 7^^,J2)^ Finally we want to make the re- 
mark that there are no interference terms between 
the Bethe-Heitler and the Compton reactions in 
(10) if one integrates over all final state momenta. 

The complexity of the second order heavy 
quark coefficient functions prohibits publishing 
them in an analytic form, except for Lf^'^'^\ 
which is given in Appendix A of ||], so that they 
are only available in large computer programs 
|l[, involving two-dimensional integrations. To 
shorten the long running time needed for the com- 
putation of the structure functions in (5) we have 
tabulated the coefficient functions in the form of a 
two dimensional array in the variables i] and ^ in 
a different computer program |^ . These variables 
are defined by 



(l-z). 



-,2 ■ 



(11) 



However when ^ 3> 1 {Q^ 3> m?) numerical in- 
stabilities appear so that it is desirable to have 
analytic expressions for the heavy quark coeffi- 
cient functions in that region. Moreover it turns 
out that for ^ > 10 the asymptotic expressions for 

(9) ^'^'^ ^2^'^"^^ (10) approach the exact ones 
so that the former can be used for charm produc- 
tion at the HERA coUidier provided > 22.5 



(k) 



X In' — , 

(fc) 



(GeV/c^) {rric = 1.5 (GeV/c). Furthermore one 
can use these asymptotic formulae in the context 
of the variable flavour number scheme as has been 
explained in [Q. 

2. METHOD 

We will now explain how the asymptotic form 
{Q^ ^ m?) of the heavy quark coefhcient func- 
tions Hi J , Li J (5) can be derived using the 
renormalization group and the operator product 
expansion (OPE) techniques. Using these tech- 
niques one can avoid the computation of the cum- 
bersome Feynman integrals and phase space in- 
tegrals which arise in the calculation of the pro- 
cesses in (8) - (10). In the limit ^ rri^ the 
heavy quark coefRcient functions behave logarith- 
mically as 

Q'^ m'^\ ^ {k).(i)( 

(12) 

with a similar expression for Lj'^''. The above 
large logarithms originate from collinear diver- 
gences which arise when is kept fixed and 
m? 0. As has been shown in Q each fixed 
order term in expression (12) can be written as 

where the power of as has to match on the left 
and right hand sides. There is a similar expres- 
sion for [i = 2, L; k,j = q,g)- Notice that 
we have suppressed the dependence on the scal- 
ing variable z in (12) and the convolution symbol 
(8> is defined by 

dz2 5{z - Z1Z2) 

x/(^i)5(^2) . (14) 
The light quark and gluon coefficient functions 
Ci^k have been calculated up to order in [^. 
The operator matrix elements (OME's) Akj are 
now also known up to the same order in pertur- 
bation theory (see [^). They are defined by 

(Til \ 
-j^) =< ]\0k\3 > , (15) 
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where Ok represent the local operators which 
show up in the operator product expansion of 
the two electromagnetic currents which appear in 
the calculation of the process (1). Notice that 
the OME's in (15) are finite which means that 
all renormalizations and mass factorizations have 
already been carried out. The last operation is 
needed because of the collinear divergences which 
appear in the OME's when the external on-mass- 
shell massless quarks and gluons are coupled to 
internal massless quanta. The ultraviolet and 
collinear divergences are regulated by using the 
method of n-dimensional regularization. The re- 
moval of these divergences has been done in the 
MS-scheme. For the computation of the heavy 
quark coefRcient functions H^^^ (8) H-'^J (9) we 

need the OME's Aq^ and respectively, which 
are given by the Feynman graphs in figs. 1,2 of 
The Bethe-Heitler coefficient functions -ffff'*'^'' 
(10) requires the calculation of ^g^'^^'' whereas 
for the Compton coefRcient function -^q^'*'^'' (10) 

one has to compute A^q'^^"^ The results for these 
OME's can be found in appendix C of Substi- 
tution of Akj and Ci^k in (13) leads to the asymp- 
totic heavy quark coefficient functions which are 
presented in Appendix D of 

3. RESULTS 

We are now interested to find out at which val- 
ues of ^ (11) or the asymptotic coefficient 
functions approach the exact ones computed in 
1^ and HJ. For that purpose we define the ratio 
Rf] which is given by 



R 



asymp,(£) 



(z,C,m2/^2) 



where Hf^^"^^ and ff^^"^^ stand for the exact 
10], ^ and asymptotic j|] heavy quark coefR- 
cient functions respectively. Choosing /i^ = 
and the range 5 < 1^ < 10^, we have plotted as 

(2) (2) 

an example Rf^ g in fig.l and i?2 g in fig. 2 for 
z = 10~^ and z = lO""*. The reason that we 
have chosen these two ratios is that the coeffi- 
cient functions Hj^ ^ and H2 j (9) constitute the 
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Figure 1. R}^'g plotted as a function of ^ for fixed 
z = 10~^ (solid line) and for z = 10~^ (dashed 
line). 



(2) 

Figure 2. ' plotted as a function of ^ for fixed 



z = 10 (solid line) and for z = 10 (dashed 
line) . 



bulk of the radiative corrections to the Born reac- 
tion (8). From fig.l we infer that and 
jj'^yj^P''^^^ coincide when ^ > ^min 
there is essentially no difference between the ra- 
tios for z = 10~^ and z = 10~^. In the case 



10^ and 



of H'^Y''^'^^ and H™^''> (see fig.2) the above 
^- value is much smaller and both coefhcient func- 
tions coincide when ^ > ^^in — 10, which is quite 
insensitive to the values chosen for z. The rea- 

(2) 

son why the convergence of R)^ ^ to one is so slow 
in comparison to i?2 ^ is unclear. Apparently the 
logarithmic terms in ^2™"=*' (2) ^^^^^ dominate 
the coefficient functions at much larger values of ^ 
than is the case for ^ similar observa- 

tion has been made for Hf^'^^K The small value 
found for ^min in the case of i?2,g is very inter- 
esting for charm production where F2(x,Q^,m\) 
can be measured with much higher accuracy than 
Fl{x,Q'^ ,rn?). Since H^^^ dominates the radia- 
tive corrections to F2{x,Q^ ,ml) one can state 
that for Q2 > 22.5 (GcV/c)^ (m^ = 1.5 GeV/c) 



rasymp,(2) 



the exact coefficient functions can be replaced by 
their asymptotic ones. However before one can 
draw definite conclusions about the dominance of 
the terms ln\Q'^ /m?) on the level of the struc- 
ture functions ones first to convolute the heavy 
quark coefficient functions with the parton densi- 
ties (see (5) ). This will be done in future work. If 
it turns out that the above logarithms also dom- 
inate Fk{x, , w^), in particular for k — 2, then 
these terms have to be resummed using the renor- 
malization group equations. This is done using 
the variable flavour number scheme approach Q . 
One of the features of this method is that one has 
to define a charm density in the proton which is 
a convolution of the OME's Akj (15) and the 
light parton densities S and G in (5). Hence for 
>> the charm quark behaves like a light 
parton provided the large logarithmic terms dom- 
inate the proton structure functions in (5). 
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